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ABSTRACT 

The MACHO and OGLE collaborations have argued that the three longest-duration bulge microlensing 
events are likely caused by nearby black holes, given the small velocities measured with microlensing 
parallax and non-detection of the lenses. However, these events may be due to lensing by more numerous 
lower-mass stars at greater distances. We find a-posteriori probabilities of 76%, 16%, and 4% that the 
three longest events are black holes, assuming a Salpeter IMF and 40 Mq cutoff for neutron-star- 
progenitors; the numbers depend strongly on the assumed mass function, but favor a black hole for 
the longest event for most standard IMFs. The longest events (> 600 days) have an a-priori ^26% 
probability to be black holes for a standard mass function. We propose a new technique for measuring 
the lens mass function using the mass distribution of long events measured with ACS, VLTI, SIM, or 
GAIA. 

Subject headings: black hole physics — gravitational lensing — Galaxy: stellar content 



1. INTRODUCTION 

Counting the number of black-hole stellar-remnants in 
our Galaxy is complicated by their faintness. Only a few 
dozen of an estimated 10 7 -10 9 black-holes are observed as 
X-ray binaries (van den Heuvel 1992). Black- hole numbers 
estimated from stellar-population synthesis are inaccurate 
due to variations in IMF with metallicity, uncertainty in 
the star-formation history, and uncertainty in the cutoff 
mass for neutron-star progenitors. To count the bulk of 
the black holes (isolated or in wide binaries) requires a 
different observable signature, the best candidate being 
gravitational microlensing since it only relies on gravity. 

Microlensing events are primarily characterized by the 
event timescale, t = 2Re/v, where Re is the Einstein ra- 
dius in the lens plane, Re = [AGMc~ 2 Dsxy] 1 / 2 , v is the 
velocity of the lens perpendicular to the observer-source 
axis, M is the lensing mass, Dl,s are the distances to the 
lens or source, x — Dl/Ds is the fraction of distance of the 
lens to the source, and y = l—x. The motion of the Earth 
causes magnification fluctuations during long events, an 
effect called "microlensing parallax" (/x-tt), allowing one 
to measure the reduced velocity, v = (vl — xvs )/y — v Q 
(Gould 1992), where L,S,Q stand for the lens, source, 
and Sun and the velocities are perpendicular to the lens- 
ing axis. Given D$ and v, the lens mass is then solely a 
function of x, 



M(x) = (vic) 2 y/(16GxD s ). 



(1) 



The mass is for x = 1, rising to oo for x = 0, thus to 
estimate M requires knowing x. 

Bennett ct al. (2001, hereafter B01) and Mao et 
al. (2002, hereafter M02) claim that the three longest- 
duration microlensing events discovered toward the bulge 
with fi — 7r measurements are likely black holes. To arrive 
at this conclusion, they assume that each lens is a member 
of a population that has a velocity and spatial distribution 
characteristic of the disk and bulge. Here we in addition 
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assume that the lenses are drawn from a population that 
has a mass function characteristic of stellar and stellar 
remnant populations in the disk and bulge, and is inde- 
pendent of x. This latter assumption is not true for young 
bright disk stars (which do not contain much of the disk 
mass), but may be true for lower mass stars and compact 
remnants in the disk and bulge. We then find that the 
probability that these events are black holes is somewhat 
reduced. 

In §2, we apply the analysis of B01 to MACHO-99-BLG- 
22 and estimate the implied total black hole number. In 
§3, we include a mass function in our prior for the lens 
mass. In §4, we consider the timescale distribution for 
bulge events and describe how mass measurements of long 
events may provide an estimate of the lens mass function. 

2. MASS ESTIMATE FROM REDUCED VELOCITY ONLY 

B01 and M02 present seven microlensing events toward 
the bulge of the Galaxy with fx— 7T measurements as shown 
in Figure 1 . These events have greatly decreased \ 2 when 
the /x-7r effect is included and are unlikely to be caused by a 
binary lens or source. B01 have computed a distance prob- 
ability for MACHO-96-BLG-5 and MACHO-98-BLG-6 us- 
ing a likelihood function which depends on the observed 
reduced velocity and assumes the source is in the bulge. 
Using equation (1) they find that the most likely masses for 
the their two longest events are > 3M Q at 68% confidence. 
Given that these lenses are not detected as main-sequence 
stars, they argue that they must be black holes. Using the 
likelihood function of B01, we find that the mass proba- 
bility for MACHO-99-BLG-22 (OGLE-1999-BUL-32) has 
M = llig 2 M© and a 81% probability of being a black 
hole using the best-fit parameters of Bennett et al. (2002, 
B02). This differs from the results of B02 due to different 
assumed disk and bulge velocity dispersions. 

If all three events are due to disk black-holes, we can 
estimate the total number. The detection efficiency, e, for 
events from 100-400 days is 20% (Alcock et al. 2000), 
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which we assume also holds at longer t (these events are 
drawn from the alert sample which have a different se- 
lection criterion, so a different e may apply). We have 
computed the timescale probability for black-hole lensing 
events, and find that only ~ 40% have t > 1 year, so £ ~ 
13%. These three events yield r ~ 5 x 10 _7 (e/0.1)" 1 , im- 
plying 7V total - 5 x 1O 8 ((M)/9M )- I (e/O.l)- 1 disk black 
holes. This estimate is larger than estimates based upon 
the expected ratio of black-hole to neutron-star remnants, 
~ 10 8 , (Shapiro & Teukolsky 1983, van den Heuvel 1992) 
and chemical enrichment by supernovae within the Milky 
Way, ~ 2 x 10 s (Samland 1998). This discrepancy indi- 
cates that either e = 20% is too low, the IMF was much 
more top-heavy in the past, or the events are due to more 
distant, low-mass stars. We next explore the third option. 
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Fig. 1 - The velocity vectors, v, and event durations, 
i (proportional to circle size), for seven /z-7r microlensing 
events from B01 and M02. The shaded regions show the 
MACHO bulge fields, while the dashed line shows the re- 
gion where the bulge density is half of the central bulge 
density in the sky plane at the distance of the bulge. Bold- 
face indicates three black hole candidate events. 

3. MASS FUNCTION PRIOR 

To recap, the small reduced velocities may indicate 
that the three longest events most likely lie within 3 kpc 
and thus should have masses > 3M . However, if a 
lens is more distant, it is less massive (equation 1) and 
thus drawn from a more abundant population, which can 
compensate for the small probability of distant lenses to 
have a small reduced velocity. Here we investigate the 
mass likelihood given v and t using the spatial probabil- 
ity dn L /dM = p(x)(j)(M)/(M), where <p(M) is the mass 
function of a 10 Gyr-old stellar population (independent 
of a;), / <j>(M)dM = 1, and (M) = J dMM(j>(M) is the 
mean stellar mass. 

To compute </>(M), we start with the broken power-law 
IMF from Kroupa (2002), and, following Gould (2000), 
convert stars with masses M cms < M < M cw d into white 
dwarfs, M cw d < M < M cns into neutron stars, and 
M cns < M into black holes. We set M cms = 1M Q and 



M cw d — 11M Q (Samland 1998), while we vary M cns . 
We use Gaussian distributions for the white-dwarf and 
neutron-star mass functions with M w d = 0.5 ± 0.1M Q 
(Bragaglia et al. 1995) and M ns = 1.35 ± 0.04M Q 
(Thorsett & Chakrabarty 1999). Black holes we describe 
with 4>bh oc M~ ' 5 for 3M < M < 15M , consistent 
with the measured mass function of black holes in X-ray 
binaries (which may have a different mass function than 
isolated black holes). We compute the relative numbers 
of each compact remnant by varying the high-mass IMF 
slope 2 < (3 < 3, where 4>imf oc M~@ , equating the num- 
ber of each compact remnant with the number of stars in 
the IMF with appropriate mass range. We assume that 
this mass function applies in both the disk and bulge. 

The number of black holes strongly depends on the slope 
of the IMF, /?, and the maximum neutron-star-progenitor 
mass, M cns . We choose as a fiducial value the Salpeter 
(3 = 2.35, consistent with the observed IMF (Kroupa 
2002), and M cns = 40M Q , consistent with chemical- 
evolution models of the Galaxy (Samland 1998). The re- 
sulting mass function has mass fractions of 7% in brown 
dwarfs, 77% in main-sequence stars, 13% in white dwarfs, 
1% in neutron stars, and 1.5% in black holes. This corre- 
sponds to a total of about 10 9 neutron stars and 2 x 10 8 
black holes in the Milky Way. 

Using Bayes' theorem, the likelihood of a lens to be 
at a given distance, x, for fixed source distance D$ is 
£(x|i,v) = £(i,v|x) for no prior on the distance and no 
errors on the timescale and velocity measurements. For a 
survey of duration T covering a solid angle of SI, then the 
solid angle cross section for magnification by > 30% by a 
lens with reduced velocity v is 2ReVI)T / D 2 L . We must then 
multiply by the total number of lenses at that distance and 
integrate over the phase space of the velocity distribution 
of lenses and sources. We include lenses and sources in the 
disk and bulge using the bar model described in Dwek et al. 
(1995) and Han & Gould (1995) and a double-exponential 
disk with scale length 3 kpc, scale height 325 pc, Galac- 
tic center distance 8 kpc, solar circular velocity 200 km/s, 
velocity ellipsoid a z : : a r = 1 : 1.3 : 2, a z = 17 km/s 
at the solar circle with a scale length of 6 kpc, and asym- 
metric drift of of/120 km/s (Buchalter, Kamionkowski & 
Rich 1997). The likelihood function becomes 



C(i,v\x) = N s / dMdv s dviJlR E vyT 



dn L 
~dM 

xD s f(vL)f(v s )5 2 (v(v s ,vi,x)-v)5(i(M,v,x)-i).(2) 

The first delta function picks out the particular v while 
the second delta function picks out the particular i. We 
convert the velocity delta function in v to a delta func- 
tion in vl , then integrate over vl and vg , and convert the 
second delta function to a delta function in mass: 



C(t, v\x) 



N s AR E D s y z Tv M(f>{M) 

i WT 



PL- 



2-K(Jl(Jb 



-,(3) 



where v = {v u v b ) = v L - xw s - y(v Q + v), af b = 

a L lb + x<la %ibi an d h b denote the components directed 
in galactic longitude and latitude, respectively. Compar- 
ing to the expression in B01, we see that the likelihood 
distributions agree if 4>{M) oc M~5. We convert from the 
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likelihood of a lens lying at a distance x to a likelihood 
in M using equation (1). We then integrate over sources 
at distance Ds using a model for the bulge and disk with 
luminosity functions measured with HST. The final con- 
straint is to require that the main-sequence lens stars not 
exceed the flux limits of B01 and M02 (we ignore the con- 
tribution from red-giant lenses). 

Figure 2 shows the computed mass probability for the 
above mass function. The white-dwarf, neutron-star, 
main-sequence cutoff, and black-hole mass cutoff show up 
as peaks in the probability distribution. The greater num- 
ber of lenses at small mass compensates for the decreased 
probability for low v events, biasing the probability to- 
ward smaller mass. The probability that each lens is a 
black hole, defined as M > 3M Q , is changed from above: 
4% for MACHO-98-BLG-6, 16% for MACHO-96-BLG- 
5, and 76% for MACHO-99-BUL-22. MACHO-98-BLG- 
6 and MACHO-96-BLG-5 are most likely main-sequence 
or white-dwarf stars at > 4 kpc; better flux limits may 
rule out a main-sequence star. These probabilities change 
by one order of magnitude depending on [3 and M cns , as 
shown in Table 1 (the labels are MACHO-XX-BLG-XX, 
while the first event is also OGLE-1999-BUL-32). The 
first two columns describe the assumed mass function, 
columns 3-5 give the black-hole probability for each event 
labeled by their MACHO event number, while the last two 
columns give the average of M 1 / 2 and M 2 for black holes 
divided by the average for all stars given the assumed IMF. 

4. PROBABILITY OF LONG-TIMESCALE LENSING 

We next compute the expected distribution of events as 
a function of timescale. For long timescales the a-priori 
differential probability distribution scales as t~ A (Mao & 
Paczynski 1996). For solar-mass stars in the disk and bulge 
this behavior occurs for t > 200 days. We can rescale the 
probability distribution for different masses 



d 2 p 
dMdt 



= (M)-^(M) 



dp(M Q ) 



(4) 



where M is measured in units of M and dp(M Q )/di is the 
timescale probability distribution assuming all lenses have 
mass Mq. Since at long timescales dp{Mo)/dt scales as 
i~ 4 , we can rescale this equation and integrate over mass 



dp 
di 



pM 2 
JM 1 



d 2 p 

dMdt 



dM oc T 



7 

J Mi 



dM<p(M)M 2 , (5) 



1/2 

for t > M 2 200 days. The fraction of all microlensing 
events as a function of mass scales as M 1 / 2 <j)(M), while 
the fraction of events in the long-timescale tail scales as 
M 2 (j>(M), strongly favoring black holes. Table 1 shows the 
average of M 1 / 2 and M 2 for black holes divided by mass- 
function average, showing that the fraction of black-hole 
events in the long-timescale tail is increased by one to two 
orders of magnitude as compared to the total number of 
events. Thus, with a measurement of the mass distribu- 
tion of events at long timescales, one can directly infer the 
mass function if it is independent of the location in phase 
space. 

Table 1: Black Hole Probabilities 
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We have computed the probability distribution for 
lensing by solar-mass stars (see Han & Gould 1995, 
1996 for the computation technique), and then con- 
volved with our assumed mass function to give the 
timescale distribution. This is compared to the ob- 
served distribution of MACHO events in Figure 3. 
We have binned the 252 alert events from MACHO 
(http://darkstar.astro.washington.edu) and then com- 
puted i 2 dn/di as Nit 2 /(ejAij), where i labels the num- 
ber of the bin, errors are Poisson. Since the efficiencies, 
ej, have not been measured for the alert events, we have 
used the efficiency curve from Alcock et al. (2000), ex- 
trapolating to longer timescale. We note that the alert 
events timescales may be affected by blending and that 
a full analysis of efficiencies is required. The rise at long 
timescale may indicate that 20% underestimates the effi- 
ciency or that the longest two events are a statistical fluke. 
An efficiency of 100% is indicated for the the last three 
bins by an arrow. Decreasing /3 to 2 improves agreement 
for the shorter timescales, but the longest event still lies 
well above the predicted value, so it may be an improbable 
event. 

We have tried to explain the rise at long timescale 
by varying the assumptions in the Galactic kinematics 
and density, but all reasonable modifications fail. To ex- 
plain this discrepancy with a different mass function re- 
quires stars of mass ~ IOOMq since the timescale scales 
as M 1 / 2 and the probability distribution peaks at ~ 100 
days for 1M Q lenses. However, 100M© lenses are ruled 
out by the /i—tt observations which indicate that the three 
longest timescale events would have to be within a few 
hundred parsecs. This would require an even larger black 
hole mass density to explain the number of these events, 
~ 1M Q pc~ 3 , which exceeds the Oort limit. Also, for 
events near the peak of the timescale distribution for 
1OOM lenses, v should be much larger than observed. 
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Fig. 3 - The distribution of event timescales, Pdn/di. 
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Fig. 2 - The differential (solid) and integrated (dashed) likelihood of M given the observed v and i. (a) MACHO-99-BLG- 
22 (OGLE-1999-BUL-32) (b) MACHO-96-BLG-5 (c) MACHO-98-BLG-6. The mass function parameters are /? = 2.35 
and M cns = 40M Q . The distance scale at the top of the plots shows Dl(M) if Ds = 8 kpc. The angle ip is the velocity 
direction in Galactic coordinates measured from the rotation direction towards the north Galactic pole. The t^-band flux 
limits are from (a) M02 (b) Bennett et al. (2002) and (c) B01. 



The histogram shows the MACHO alert events, while the 
lines show the predicted distribution for all masses (solid) , 
neutron stars (dashed), and black holes (dotted) for the 
fiducial IMF. 

5. CONCLUSIONS 

We have shown that including the timescale as a con- 
straint on the black hole probability may change the con- 
clusion as to whether or not the lens is a black hole. The 
probability is most robust for MACHO-99-BLG-22, which 
is most likely a black hole for the range of mass functions 
we have explored, while the probability for MACHO-98- 
BLG-6 and MACHO-96-BLG-5 is strongly dependent on 
the assumed mass function. A fi-ir measurement does not 
provide sufficient information to estimate the mass of a 
given event, but may result in interesting limits on the 
mass. 

The contribution of a given population of stars to the 
long-timescale tail (£ > 200M 1 ' 2 days) of the lensing dis- 
tribution depends on (M 2 ) for that population. Due to 
their larger mass, ~ 26% of events with t > 600 days 
should be black-hole lenses for j3 — 2.35 and M cns = 
40M Q ; this fraction may be reduced if black holes have 
a velocity dispersion much larger than that of disk and 
bulge stars. To constrain (3 and M cns with the timescale 
distribution will require more events and a careful esti- 
mate of the microlens detection efficiency as a function of 
t. The parallax events may have an increased efficiency 
since they are less affected by observing gaps and have a 
wider cross section due to the Earth's motion (B01). 

The OGLE III project may detect 10 3 microlensing 
events toward the bulge each year (Paczyhski, priv. 
comm.), which should include ~ 3 black- hole lensing 
events and ~ 5 neutron-star events for (3 = 2.35 and 
M cns — 40M Q . Since the timescale of black-hole events 



is longer, ji— tt measurements are feasible from the ground 
(as demonstrated by B01). If these can be followed up 
with astrometric observations, then the mass can be de- 
termined when combined with /i-7r information (Paczyhski 
1998, Gould 2001, Boden, Shao, & van Buren 1998). Since 
the typical mass of black holes, ~ 9M Q , is much higher 
than the typical lens mass, ~ 0.5M Q , the typical astro- 
metric signal for black-hole lens events with a source in 
the bulge will be ~ 3 mas which will easily be measured 
with ACS (M02), VLTI (by 2004, Delplancke et al. 2001), 
SIM (by 2009), or GAIA (by 2010-12, Belokurov & Evans 
2002), compared to 0.7 mas for a typical main-sequence 
lens. The deviation of the centroid is 59 = 9eu/(u 2 + 2), 
where u is the normalized impact parameter as a function 
of time and Oe is the Einstein angle (Paczyhski 1998). 
Since u is determined by the photometric light curve, ob- 
servations of 68 at two different u will suffice to determine 
0e- Combined with Re = vi/2 measured photometrically 
from the /i-n will allow a determination of the mass of 
the lens, M = c 2 ReOe/{^G). With mass measurements 
for many long-timescale events, the mass function can be 
deduced since the probability for lensing is proportional 
to i- 4 M 2 (j>(M) for a phase -space distribution that is inde- 
pendent of mass (equation 3). The M 2 dependence means 
that detection of black holes will be favored so that one 
can constrain f3 and M cns . 

Finally, direct observations of the lenses provide another 
avenue toward determining their nature. If they are black 
holes, they may accrete from the interstellar medium or 
from a companion wind. With plausible assumptions, the 
X-ray radiation could be detectable with current space- 
based observatories (Agol & Kamionkowski 2001). 
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